Probing the intrinsic Josephson potential in Bi-2212 by thermal activation 
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We study thermal fluctuation phenomena in small Bi-2212 intrinsic Josephson junctions. Being 
able to measure switching currents of a single intrinsic junction, we observe that it's statistics can 
be very well described by thermal activation from a periodic Josephson potential with the sinusoidal 
current-phase relation. This is a direct evidence for the dc-intrinsic Josephson effect and the first 
unambiguous confirmation of the tunnelling nature of interlayer transport in strongly anisotropic 
high temperature superconductors. Furthermore, the fluctuation-free critical current, extracted 
from the analysis of switching current statistics, exhibits a temperature dependence typical for 
superconductor- insulator- superconductor tunnel junctions. 
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The nature of interlayer transport in high super- 
conductors (HTSC) has been a long standing question 
P|. It is established that in extreme anisotropic Bi- and 
Tl-based HTSC intrinsic Josephson effect exists between 
superconducting Cu02 planes 2]. Both dc and ac intrin- 
sic Josephson effects were observed in Bi2Sr2CaCu208+a; 
(Bi-2212), in the form of flux quantization Fiske 
|3| and Shapiro 6] steps in Current- Voltage character- 
istics (IVC's) and the Josephson plasma resonance 0- 
On the other hand, the role of blocking Bi-laycrs and the 
nature of interlayer coupling in Bi-2212 is still unclear, 
partly due to difficulties in analysing intrinsic Joseph- 
son data caused by stacking and strong electromagnetic 
coupling of intrinsic Josephson junctions (IJJ's) j^j- 

The type of Josephson coupling can be deduced from 
the dependence of the Josephson energy Ej on the phase 
difference tp, or, similarly, from the Josephson current- 
phase relationship Is{<p) = (2e//i)9i?j/9(^, which varies 
from a saw-tooth like for metallic weak links to the si- 
nusoidal for superconductor- insulator- superconductor 
(SIS) tunnel junctions 0. The Ej{ip) determines the 
junction electrodynamics, which is equivalent to motion 
of a particle in a "tilted wash-board" potential, cre- 
ated by superposition of the periodic Josephson poten- 
tial Ej{ip) and the work done by the current source, 
— {h/2e)I(p. At a finite temperature, T, the particle can 
escape from the potential well as a result of thermal fluc- 
tuations, see Fig. 3 a). This corresponds to switching 
of the junction from the superconducting to the resistive 
state. The rate of thermal escape 
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is a sensitive probe of the Ej{(p) both via the attempt fre- 
quency, LUa, i.e., the frequency of oscillations at the bot- 
tom of the potential well, and, particularly, via strong ex- 
ponential dependence of Ft on the potential barrier AU. 
Thus, the switching current statistics carries direct infor- 



mation about the shape of Ej{tp) and, therefore, about 
the nature of Josephson coupling. 

In this letter we study the effect of thermal fluctua- 
tions in small Bi-2212 intrinsic Josephson junctions. Be- 
ing able to measure switching currents of a single IJJ in 
an electrically shielded environment, we observe that it's 
statistics can be very well, and without fitting param- 
eters, described by thermal activation from the tilted 
wash-board potential with the sinusoidal current-phase 
relation. This is a direct evidence for the dc-intrinsic 
Josephson effect and the first unambiguous confirmation 
of the tunnelling nature of interlayer transport in Bi- 
2212. We also demonstrate that thermal fluctuations 
dramatically affect properties of small IJJ's, resulting in 
strong suppression of the switching current and unusual 
temperature dependence in the whole T— range. How- 
ever, the fluctuation- free critical current /co, extracted 
from the analysis of switching current histograms, ex- 
hibit T— dependence typical for SIS junctions, consistent 
with tunnelling nature of the interlayer coupling. 

The IJJ's were fabricated by etching small mesa struc- 
tures on top of Bi-2212 single crystals. We developed a 
simple procedure capable of fabricating deep sub-micron 
multi-terminal IJJ's. Fig. 1 shows a sketch of fabrica- 
tion procedure. It involved self-alignment cross-bar pho- 
tolithography, during which an insulating CaF2 layer was 
formed using lift-off, see Fig. 1 a), and ~ 3 x 3 — 5 x 5fim^ 
mesas were formed in a self-aligned manner, Fig. 1 c), 
at the crossing between the narrow-long mesa, Fig. 1 a) 
and bar- like electrodes. Fig. 1 b). Finally, the sample was 
transferred into standard Focused Ion Beam (FIB) sys- 
tem (FEI Inc. FIB-200), and a smaller mesa was trimmed 
by cutting of a portion of the mesa, as shown in Fig. 1 
d). Due to self-ahgnment at the previous stage there is 
no parasitic area below the electrode and deep-submicron 
mesas can be fabricated. Fig.l e) shows IVC's of mesas 
before and after FIB trimming to sub- micron dimensions. 
Both IVC's exhibit a knee at the sum-gap voltage, fol- 
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FIG. 1: A sketch of sample fabrication: a-c) Self - align- 
ment cross-bar lithography, d) secondary electron image of 
the sample after FIB trimming and e) I-V curves of mesas 
before and after FIB trimming. 



lowed by almost temperature independent high bias re- 
sistance [l^ . which is typical for SIS tunnel junctions. 
The increase in resistance after trimming is in agreement 
with ^ 80 fold decrease in the mesa area. 

Fig. 2 shows JVC's, normalized by the mesa area, at 
T ~ 6K for the initial 6 x Z^m? mesa and two smaller 
mesas obtained by cutting the initial mesa in two parts 
with the FIB. The multi-branch structure of the IVC's is 
due to one-by-one switching of the IJJ's from the super- 
conducting to the resistive (quasiparticle) state. Naively, 
it could be expected that the normalized IVC's should 
collapse in one, since the critical current density is a ma- 
terial property, independent of the junction area. In re- 
ality, the measured switching current density decreases 
with the mesa area. To understand whether this is caused 
by deterioration of the IJJ's during FIB trimming and 
thermal cycling, we checked the scaling of quasi-particle 
branches in the IVC's. The thin solid lines in Fig. 2 rep- 
resent multiple-integers of the fit to the first branch of the 
initial mesa. It is seen that quasiparticle branches of the 
initial mesa are perfectly periodic. However, it is also 
seen that the normalized quasiparticle branches of the 
trimmed mesas fall onto the same fitting curves. This im- 
plies that there is no visible deterioration of IJJ's during 
FIB trimming and thermal cycling. Therefore, the ob- 
served decrease of the switching current density in small 
mesas should be attributed to enhanced thermal activa- 
tions due to decrease of the potential barrier A[/, which 
is proportional to the Josephson energy _Bjo = {hc/2e)IcQ 
and scales with the junction area. 

For a quantitative analysis of the switching currents in 
small mesas we consider the effect of thermal activation 
from a tilted wash board potential. The probability to 



FIG. 2: Current density vs. Voltage curves for a mesa before 
and after FIB cutting in two non-equal parts. The measured 
switching current density decreases with the size of the mesa. 
Solid lines are fits to quasiparticle branches, which are mul- 
tiple integers of the first branch. It is seen that quasiparticle 
branches of all junctions fall onto the same line irrespective 
of the size of the mesa, which implies that there is no junc- 
tion deterioration. Horizontal dashed lines show values of the 
fluctuation free critical current density (the top line) and the 
most probable switching current densities for the three mesas, 
obtained from the fit shown in the inset. Inset: The solid line 
shows simulated dependence of the most probable switching 
current Imax vs. T/Ejo, calculated for a classical thermal 
escape from a tilted wash board potential. Circles indicate 
the best fit for the there mesas. 



measure the switching current / is given by 
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where dl/dt is the current sweeping rate. For under- 
damped junctions, Q = LUpRC 3> 1, the prefactor in 
Eq.lHJ is at = A/ [{I + QkBT/l.8AU)°-^ + l]"^ [ij. In 
simulations we will assume the sinusoidal Is{f), for which 
M^a= Up{l-{I/IcO?f'\ where = {2eho/hcCf/^ 
is the Josephson plasma frequency, and 
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The inset in Fig. 2 shows the simulated dependence of 
the most probable switching current Imax as a function of 
the ratio T/Ejo, made for the typical experimental con- 
ditions. It is seen that Imax decreases with Ejo, which 
in turn is proportional to the area of the junction. The 
circles in the inset show the fit to this universal depen- 
dence made using three fitting variables: the fluctuation 
free critical current density Jco (the same for all mesas), 
the effective noise temperature T'e// for the initial mesa 
and Teff for trimmed mesas (the same for both since 
they were measured in the same run). The horizontal 



3 




30 40 50 60 70 80 90 100 110 120 



/(^lA) 




-1 1 , 2 3 20 40 60 80 

9 / 7C J (K) 

FIG. 3: a) The shape of tilted wash-board potential for 
I /IcO ~ 0.5. b) Switching current histograms of a single IJJ 
at different temperatures. Solid lines represent fits to classi- 
cal thermal escape from the tilted wash-board potential. The 
inset shows the quality of the fit for T — 46.2K. Dots and 
lines represent experimental data and fits, respectively, c) 
The effective "escape temperature" vs. temperature. 

dashed lines in Fig. 2 show J^o (the top hne) and the 
most probable switching current densities for the three 
mesas, obtained from such a fit. Good agreement is seen 
between fitted and measured switching currents for all 
three mesas. Taking into account that there was no free 
fitting parameters (three points were fitted with three 
variables) and that Te// = lOK and 12K are similar 
and only few degrees above the substrate temperature 
T ~ 6K, we conclude that the observed decrease of the 
switching current density in smaller mesas is the result 
of enhanced thermal fluctuation caused by reduction of 
the Josephson coupling energy Ejq with junction area. 

The thermal escape rate, Eq.(l), strongly depends on 
the shape of the wash board potential, AU{ip), which in 
turn depends on the current-phase relation Js{^)- There- 
fore, the probability distribution of the switching current 
P{I) contains explicit information about the shape of 
Josephson potential Ej{ip). Since Fig. 2 b) indicates 
that the switching current density of small IJJ's is well 
described by thermal activation from the wash-board po- 
tential, we should be able to probe the intrinsic Joseph- 
son potential by studying switching current statistics. 
Unfortunately, previous studies of the switching current 
statistics in Bi-2212 mesas revealed that switching cur- 
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FIG. 4; Temperature dependencies of the most probable 
switching current Imax (circles), the retrapping current Ir 
(rhombs), the extracted fluctuation-free critical current Ico 
(squares), and the superconducting gap A (triangles) for the 
same IJJ as in Fig. 3. It is seen that Imax has unusual lin- 
ear dependence in the whole T— range. However, Ico{T) is 
normal, close to A(T) and consistent with tunnelling nature 
of the interlayer transport. For comparison, T— dependencies 
of the energy gap A{BCS) and the critical current Ic{AB) 
for conventional SIS junctions are shown by solid and dashed 
lines, respectively. 

rent histograms of mesas containing several stacked IJJ's 
may be very unusual. It was reported that histograms 
of stacked IJJ's may contain multiple peaks £^nd 
appeared to be extremely broad |3, , up to ~ten 

times broader than expected. Such anomalous behav- 
ior was attributed to the presence of multiple metastable 
states (fluxon modes), which appear due to coupling of 
junctions in the stack The existence of metastable 

states results in a multiple valued critical current and 
dramatic enhancement of thermal fluctuations in stacked 
Josephson junctions. Note that such behavior is not spe- 
cific to Bi-2212 mesas, but was also observed for low-Tc 
stacked Josephson junctions j^. 

In order to study the intrinsic Josephson potential we 
have to avoid metastable states. To solve this problem 
here we have studied switching of a single IJJ. The stable 
switching of a single junction can be obtained when there 
is a spread in critical currents between IJJ's in the mesa. 
Such a spread is often observed in mesas obtained by wet 
chemical etching, see Fig. 2, and is most probably caused 
by variation of the junction area due to undercut. 

Fig. 3 represents the switching current statistics of 
a single IJJ in an optimally doped Bi-2212 mesa. This 
junction had ^ 20% smaller critical current than the rest 
of IJJ's in the mesa, which was sufficient for achieving 
stable biasing without switching the rest of IJJ's. Mea- 
surements were done in a shielded room environment us- 
ing a sample-and-hold setup with the effective noise tem- 
perature ~ lOOmK T^. Fig. 3 b) shows switching cur- 
rent histograms obtained from 10240 switching events at 
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different T. The solid lines represent fits to classical ther- 
mal escape, Eqs. (1-3). The fit was made using follow- 
ing parameters: the experimental sweeping rate dl/dt — 
24:.3mA/s, the specific capacitance C = 
and the junction resistance, R, extracted from the high 
bias resistivity pc = 25flcm. Since the probability dis- 
tribution is only slightly dependent on C and R, those 
parameters were fixed during the fit to avoid ambiguity. 
The only remaining fitting parameter was the effective 
"escape" temperature T^sc, which in the absence of noise 
should coincide with T. The values of T^sc obtained from 
the fit are plotted as a function of T in Fig. 3 c). It is 
seen that for T < 72K, Tesc follows T. This, confirms 
the validity of the fitting procedure and clearly demon- 
strates the sinusoidal current-phase relation and the cos- 
inusoidal dependence of the intrinsic Josephson potential 
Ej(ip). The remarkable accuracy with which the sinu- 
soidal current-phase relationship is satisfied can be seen 
from the excellent quality of the fit, which is shown in 
detail in inset to Fig. 3 b). Such behavior is character- 
istic for high quality SIS tunnel junctions. This is the 
first unambiguous evidence for the tunnelling nature of 
interlayer transport in Bi-2212. 

At high temperatures, T > 75K, T^sc starts to de- 
crease and eventually vanishes close to Tc ~ QZK. The 
surprising collapse of thermal fluctuations close to Tc is 
associated with the change in the shape of switching his- 
tograms, which lose their characteristic asymmetric form 
and become symmetric and narrow. It is also associ- 
ated with the decrease and collapse at T ~ ^QK of the 
hysteresis in IVC, as seen from comparison of the most 
probable switching current, Imax, and the retrapping cur- 
rent /r, shown in Fig. 4. We note that the switching 
current remains sharply defined and there is no indica- 
tion for a phase-diffusion 0| up to ^ QQK. We be- 
lieve that the collapse of T^sc is caused by entering high 
dissipation regime close to Tc, in which the spread in 
switching currents is reduced by enhanced probability of 
retrapping of the rolling particle in the wash board po- 
tential. We emphasize that such unusual behavior is not 
unique for Bi-2212 IJJ's but was also observed in low- 
Tc superconductor- normal metal- superconductor junc- 
tions. Therefore, this phenomenon is not essential for the 
present work and will be discussed elsewhere 16J. 

Fig. 4 shows temperature dependence of the most 
probable switching current Imax (circles), the retrap- 
ping current 1^ (rhombs), and /co, obtained from fit- 
ting switching current histograms (squares), for the same 
IJJ as in Fig. 3. It is seen that Imax has an un- 
usual linear dependence in the whole T— range. How- 
ever, the T— dependence of Jco is quite normal and close 
to T— dependence of the superconducting gap A (trian- 
gles), which was obtained from the sum-gap knee in IVC's 
at higher bias For comparison, T— dependencies 

of the conventional ECS energy gap, /S.{BCS) and the 
Ambeokar-Baratoff value of the critical current Ic{AB) 



for conventional SIS junctions are shown in Fig. 4 by 
solid and dashed lines, respectively. It is seen that 
Ic{AB) (X A{BCS) at T < Tc/2. The experimental 
A(T) deviates somewhat from Ico{T) in the intermediate 
T— range. The deviation is likely a result of self-heating 
at the large sum-gap voltage V = 2NA/e, where N is the 
number of IJJ's in the mesa. Such deviation is in agree- 
ment with both numerical simulations |^| and in-situ 
measurement of self- heating in our mesas |l9j . Therefore, 
the unusual T— dependence of the switching current is 
solely due to thermal fluctuations, while T— dependence 
of the extracted fluctuation free critical current is consis- 
tent with the tunnelling nature of interlayer transport. 

In conclusion, having studied the switching current 
statistics of a single Bi-2212 intrinsic Josephson junction, 
we observed that it can be very well, and without fit- 
ting parameters, described by thermal activation from a 
tilted wash-board potential with the sinusoidal current- 
phase relation. This is direct evidence for the dc-intrinsic 
Josephson effect and the first unambiguous confirma- 
tion of tunnelling nature of the interlayer transport in 
strongly anisotropic HTSC. We demonstrated that ther- 
mal fluctuations dramatically affect properties of small 
IJJ's, resulting in strong suppression of the switching 
current density and unusual T— dependence in the whole 
T— range. However, fluctuation- free Jco, extracted from 
the analysis of switching current histograms, exhibit a 
T— dependence typical for SIS tunnel junctions, also con- 
firming tunnelling nature of the interlayer coupling. 
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